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Abstract. Phylogenetic networks are a generalization of phylogenetic trees that 
allow for the representation of non-treelike evolutionary events, like recombina- 
tion, hybridization, or lateral gene transfer. While much progress has been made 
to find practical algorithms for reconstructing a phylogenetic network from a set of 
sequences, all attempts to endorse a class of phylogenetic networks (strictly extend- 
ing the class of phylogenetic trees) with a well-founded distance measure have, to 
the best of our knowledge, failed so far. In this paper, we present and study a new 
meaningful class of phylogenetic networks, called tree-child phylogenetic networks, 
and we provide an injective representation of these networks as multisets of vectors 
of natural numbers, their path multiplicity vectors. We then use this representa- 
tion to define a distance on this class that extends the well-known Robinson-Foulds 
distance for phylogenetic trees, and to give an aligmnent method for pairs of net- 
works in this class. Simple, polynomial algorithms for reconstructing a tree-child 
phylogenetic network from its path multiplicity vectors, for computing the distance 
between two tree-child phylogenetic networks, and for aligning a pair of tree-child 
phylogenetic networks, are provided. They have been implemented as a Perl pack- 
age and a Java applet, and they are available at the Supplementary Material web 
page. 



1 Introduction 

Phylogenetic networks have been studied over the last years as a richer model of 

the evolutionary history of sets of organisms than phylogenetic trees, because they 
take not only mutation events but also recombination, hybridization, and lateral 
gene transfer events into account. 

The problem of reconstructing a phylogenetic network with the least possible 
number of recombination events is NP-hard [41], and much effort has been devoted 
to bounding the number of recombination events needed to explain the evolutionary 
history of a set of sequences [2, 26, 38]. On the other hand, much progress has been 
made to find practical algorithms for reconstructing a phylogenetic network from 
a set of sequences [10, 11, 23, 29, 31, 38]. 

Since different reconstruction methods applied to the same sequences, or a 
single method applied to different sequences, may yield different phylogenetic net- 
works for a given set of species, a sound measure to compare phylogenetic networks 
becomes necessary [30] . The comparison of phylogenetic networks is also needed in 
the assessment of phylogenetic reconstruction methods [21], and it will be required 
to perform queries on the future databases of phylogenetic networks [34]. 



Many metrics for the comparison of phylogenetic trees are known, including 
the Robinson- Foulds metric [36], the nearest-neighbor interchange metric [42], the 
subtree transfer distance [1], the quartet metric [9], and the metric from the nodal 
distance algorithm [6]. But, to our knowledge, only one metric (up to small varia- 
tions) for phylogenetic networks has been proposed so far. It is the so-called error, 
or tripartition, metric, developed by Moret, Nakhleh, Warnow and collaborators 
in a series of papers devoted to the study of reconstructibility of phylogenetic net- 
works [18, 19, 22, 23, 27, 28, 30], and which we recall in §2.4 below. Unfortunately, it 
turns out that, even in its strongest form [23], this error metric never distinguishes 
all pairs of phylogenetic networks that, according to its authors, are distinguish- 
able: see [7] for a discussion of the error metric's downsides. 

The main goal of this paper is to introduce a metric on a restricted, but mean- 
ingful, class of phylogenetic networks: the tree-child phylogenetic networks. These 
are the phylogenetic networks where every non-extant species has some descendant 
through mutation. This is a slightly more restricted class of phylogenetic networks 
than the tree-sibling ones (see §2.3) where one of the versions of the error met- 
ric was defined. Tree-child phylogenetic networks include galled trees [10, 11] as a 
particular case, and they have been recently proposed by S. J Wilson as the class 
where meaningful phylogenetic networks should be searched [43]. 

We prove that each tree-child phylogenetic network with n leaves can be singled 
out, up to isomorphisms, among all tree-child phylogenetic networks with n leaves 
by means of a finite multisubset of N"". This multiset of vectors consists of the 
path multiplicity vectors, or fi-vectors for short, fi(v) of all nodes v of the network: 
for every node v, ii{v) is the vector listing the number of paths from v to each 
one of he leaves of the network. We present a simple polynomial time algorithm 
for reconstructing a tree-child phylogenetic network from the knowledge of this 
multiset. 

This injective representation of tree-child phylogenetic networks as multisub- 
sets of vectors of natural numbers allows us to define a metric on any class of 
tree-child phylogenetic networks with the same leaves as simply the symmetric 
difference of the path multiplicity vectors multisets. This metric, which we call 
^1- distance, extends to tree-child phylogenetic networks the Robinson- Foulds met- 
ric for phylogenetic trees, and it satisfies the axioms of distances, including the 
separation axiom (non-isomorphic phylogenetic networks are at non-zero distance) 
and the triangle inequality. 

The properties of the path multiplicity representation of tree-child phylogenetic 
networks allow us also to define an alignment method for them. Our algorithm 
outputs an injective matching from the network with less nodes into the other 
network that minimizes in some specific sense the difference between the jU-vectors 
of the matched nodes. Although several alignment methods for phylogenetic trees 
are known [25, 32, 33], this is to our knowledge the first one that can be applied to 
a larger class of phylogenetic networks. 
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Wc have implemented our algorithms to recover a tree-child phylogenetic net- 
work from its data multiplicity representation and to compute the //-distance, to- 
gether with other related algorithms (like for instance the systematic and efficient 
generation of all tree-child phylogenetic networks with a given number of leaves), 
in a Perl package which is available at the Supplementary Material web page. We 
have also implemented our alignment method as a Java applet which can be run 
interactively at the aforementioned web page. 

The plan of the rest of the paper is as follows. In Section 2 we gather some 
preliminary material: we fix some notations and conventions on directed acyclic 
graphs, and we recall several notions related to phylogenetic trees and networks, 
the Robinson-Foulds metric for the former and the tripartition metric for the lat- 
ter. In Section 3 we introduce the tree-child phylogenetic networks and wc study 
some of their basic properties. In Section 4 we introduce the path multiplicity 
representation of networks and we prove that it singles out tree-child phylogenetic 
networks up to isomorphism. Then, in Section 5 we define and study the /x-distance 
for tree-child phylogenetic networks with the same number of leaves, and in Sec- 
tion 6 we present our alignment method. The paper ends with a short Conclusion 
section. 

2 Preliminaries 
2.1 DAGs 

Let N = {V, E) be a directed acyclic graph (DAG). We denote by di{u) and do{u) 
the in-degree and out-degree, respectively, of a node u £ V. 

A node v £V is a leaf if do{v) = 0, and internal if do{v) > 0; a root if di{v) = 0; 
a tree node if di{v) ^ 1, and a hybrid node if di{v) > 1. We denote by Vl, Vt, 
and Vh the sets of leaves, of tree nodes, and of hybrid nodes of N, respectively. A 
DAG is said to be rooted when it has only one root. 

Given an arc (u, v) G E, we call the node u its tail and the node v its head. 
An arc {u, v) G E is a tree arc if u is a tree node, and a hybridization arc if v is 
hybrid. We denote by Et and E^ the sets of tree arcs and of hybridization arcs, 
respectively. 

A node v G V is a child of li G y if {u, v) G V; we also say that u is a parent 
of V. For every node u &V, let child(ti) denote the set of its children. All children 
of the same node are said to be siblings of each other. The tree children of a node 
u are its children that are tree nodes. 

A DAG is binary when all its internal tree nodes have out-degree 2 and all its 
hybrid nodes have in-degree 2 and out-degree 1. 

Let S be any finite set of labels. We say that the DAG N is labeled in S, or that 
it is an S-DAG, for short, when its leaves are bijectively labeled by elements of S. 
Two DAGs A^, N' labeled in S are isom^orphic, in symbols N = N' , when they are 
isomorphic as directed graphs and the isomorphism preserves the leaves' labels. 
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In this paper we shall always assume, usually without any further notice, that 
the DAGs appearing in it are labeled in some set S, and we shall always identify, 
usually without any further notice either, each leaf of a DAG with its label in S. 

A path in jV is a sequence of nodes {vq, vi, . . . , Vk) such that ■Uj) G E for all 

i = 1, . . . ,k. We say that such a path starts in vq, passes through vi, . . . , v^-i and 
ends in Vk] consistently, we call vq the origin of the path, vi, . . . , Vk-i its interme- 
diate nodes, and Vk its end. The position of the node Vi in the path (vq, vi, . . . , Vk) 
is ? + 1. The length of the path {vo,vi, . . . ,Vk) is k, and it is non-trivial ii k ^ 1: 
a trivial path is, then, simply a node. We denote by u-^'y any path with origin u 
and end v. 

The height of a node is the length of a longest path starting in the node and 
ending in a leaf. 

We shall say that a path u~^v is contained in, or that it is a subpath of, a path 
u' v' when there exist paths u' ■^u and v^v' such that the path u' ^ v' is the 
concatenation of the paths u' -^u, u-^v, and v-^v'. 

A path is elementary when its origin has out-degree 1 and all its intermediate 
nodes have in and out-degree 1. 

The relation ^ on F defined by 

V there exists a path u-^v 

is a partial order, called the path ordering on N. Whenever u ^ v, we shall say 
that V is a descendant of u and also that u is an ancestor of v. For every node 
u £ V", we shall denote by C{u) the set of all its descendants, and by Cl{u) the 
set of leaves that are descendants of u: we call Cl{u) the cluster of u. 

A node u of AT is a strict descendant of a node u if it is a descendant of it, and 
every path from a root of to f contains the node u: in particular, we understand 
every node as a strict descendant of itself. For every node u & V, we shall denote 
by A{u) the set of all its strict descendants, and by Al{u) the set of leaves that 
are strict descendants of u: we call Al{u) the strict cluster of u. 

A tree path is a non-trivial path such that its end and all its intermediate nodes 
are tree nodes. A node is a tree descendant of a node u when there exists a tree 
path from u to v. For every node u eV, we shall denote by T{u) the set of all its 
tree descendants, and by Tl{u) the set of leaves that are tree descendants of u: we 
call Ti{u) the tree cluster of u. 

We recall from [7] the following two easy results, which will be used several 
times in the next sections. 

Lemma 1. Let u^v be a tree path. Then, for every other path w-^v ending in 
V, it is either contained in u-^v or it contains u-^v. 

CoroUeiry 1. If v G T{u), then v G A{u) and the path u-^v is unique. 



4 



2.2 The Robinson-Foulds metric on phylogenetic trees 

A phylogenetic tree on a set S of taxa is a rooted tree without out-degree 1 nodes 
with its leaves labeled bijectively in S, i.e., a rooted S-DAG with neither hybrid 
nodes nor out-degree 1 nodes. 

Every arc e = (u, v) of a phylogenetic tree T = (F, E) on S defines a bipartition 
oiS 

7r{e) = {Cl{v),S\Cl{v)). 
Let 7r(r) denote the set of all these bipartitions: 

7r(r) = {7r(e) \ ee E}. 

The Robinson-Foulds metric [36] between two phylogenetic trees T and T' on 
the same set S of taxa is defined as 

dflf(r,T') = |7r(T)A7r(rO|, 

where A denotes the symmetric difference of sets. 

The Robinson-Foulds metric is a true distance for phylogenetic trees, in the 
sense that it satisfies the axioms of distances up to isomorphisms: for every phylo- 
genetic trees T, T' , T" on the same set S of taxa, 

(a) Non-negativity: dfip{T,T') ^ 

(b) Separation: diiF{T, T') = if and only if T = T' 

(c) Symmetry: dRF{T,T') = dRF{T',T) 

(d) Triangle inequality: dRF{T, T') ^ dRF{T, T") + dRF{T", T') 

2.3 Phylogenetic networks 

A natural model for describing an evolutionary history is a directed acyclic graph 
(DAG for short) whose arcs represent the relation parent-child. Such a DAG will 
satisfy some specific features depending on the nature and properties of this rela- 
tion. For instance, if we assume the existence of a common ancestor of all individ- 
uals under consideration, then the DAG will be rooted: it will have only one root. 
If, moreover, the evolutionary history to be described is driven only by mutation 
events, and hence every individual has only one parent, then the DAG will be a 
tree. In this line of thought, a phylogenetic network is defined formally as a rooted 
DAG with some specific features that are suited to model evolution under muta- 
tion and recombination, but the exact definition varies from paper to paper: see, 
for instance, [3, 12-14, 19, 37, 39, 40]. 

For instance, Moret, Nakhleh, Warnow and collaborators have proposed several 
slightly different definitions of phylogenetic networks [18, 19, 22, 23, 27, 28]. To re- 
call one of them, in [18] a model phylogenetic network on a set S of taxa is defined 
as a rooted S-DAG satisfying the following conditions: 
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(1.1) The root and all internal tree nodes have out-degree 2. All hybrid nodes have 
out-degree 1 , and they can only have in-degree 2 (allo-polyploid hybrid nodes) 
or 1 (auto-polyploid hybrid nodes). 

(1.2) The child of a hybrid node is always a tree node. 

(1.3) Time consistency: If x,y are two nodes for which there exists a sequence of 
nodes {vq, vi, . . . , vj-) with vq = x and = y such that: 

— for every i = 0, . . . ,k — 1, either [vi, u^+i) is an arc of N, or {vi+i,Vi) is a 
hybridization arc of A'^, 

— at least one pair (vi, Vi-^^i) is a tree arc of 

then X and y cannot have a hybrid child in common. 

(This time compatibility condition (1.3) is equivalent to the existence of a temporal 
representation of the network [5,20]: an assignation of times to the nodes of the 
network that strictly increases on tree arcs and so that the parents of each hybrid 
node coexist in time. See [5, Thm. 3] or [7, Prop. 1] for a proof of this equivalence.) 

On the other hand, these authors define in loc. cit. a reconstructible phylo- 
genetic network as a rooted 5-DAG where the previous conditions are relaxed 
as follows: tree nodes can have any out-degree greater than 1; hybrid nodes can 
have any in-degree greater than 1 and any out-degree greater than 0; hybrid nodes 
can have hybrid children; and the time consistency need not hold any longer. So, 
reconstructible phylogenetic networks in this sense are simply rooted DAGs with 
neither out-degree 1 tree nodes nor hybrid leaves. These model and reconstructible 
phylogenetic networks are used, for instance, in [30]. 

A generalization of reconstructible phylogenetic networks are the hybrid phy- 
logenies of [4]: rooted <S-DAGs without out-degree 1 tree nodes. But although 
out-degree 1 tree nodes cannot be reconstructed, they can be useful both from the 
biological point of view, to include auto-polyploidy in the model, as well as from 
the formal point of view, to restore time compatibility and the impossibility of 
successive hybridizations in reconstructed phylogenetic networks [23, Fig. 13]. 

In papers on phylogenetic networks it is usual to impose extra assumptions to 
the structure of the network, in order to narrow the output space of reconstruc- 
tion algorithms or to guarantee certain desired properties. For instance, Nakhleh 
imposes in his PhD Thesis [27] the tree-sibling^ condition to the phylogenetic net- 
works defined above: every hybrid node must have at least one sibling that is a tree 
node. Although this condition is imposed therein to try to guarantee that the error 
metric considered in that work satisfies the separation axiom of distances (see the 
next subsection), it has also appeared under a different characterization in some 
papers devoted to phylogenetic network reconstruction algorithms [15, 16]. Indeed, 
the phylogenetic networks considered in these papers are obtained by adding hy- 
bridization arcs to a phylogenetic tree by repeating the following procedure: 

1. choose pairs of arcs {ui,vi) and {u2,V2) in the tree; 

^ Nakhleh uses the term class I to refer to these networks, but for consistency with the notations 
we introduce in the next section, we have renamed them here. 
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2. split the first into (ui,wi) and (wi,vi), with wi a new (tree) node; 

3. spht the second one into {u2,W2) and (u>2,f2), with 102 a new (hybrid) node; 

4. add a new arc {wi,W2)- 

It is not difficult to prove that the phylogenetic networks obtained in this way 
are tree-sibling, and that the binary tree-sibling phylogenetic networks are exactly 
those obtained by applying this procedure to binary phylogenetic trees. 

An even stronger condition is the one imposed on galled trees [10,11,41]: no 
tree node has out-degree 1, all hybrid nodes have in-degree 2, and no arc belongs 
to two recombination cycles. Here, by a recombination cycle we mean a pair of two 
paths with the same origin and end and no intermediate node in common. In the 
aforementioned papers these galled trees need not satisfy the time compatibility 
condition, but in other works they are imposed to satisfy it [27, 28, 31]. 

2.4 Previous work on metrics for phylogenetic networks 

While many metrics for phylogenetic trees have been introduced and implemented 
in the literature (see, for instance, [8, 35] and the references therein), to our knowl- 
edge the only similarity measures for phylogenetic networks proposed so far are due 
to Moret, Nakhleh, Warnow and collaborators in the series of papers quoted in the 
last subsection, where they are applied in the assessment of phylogenetic network 
reconstruction algorithms. We briefly recall these measures in this subsection. 

The error, or tripartition, metric is a natural generalization to networks of the 
Robinson-Foulds metric for phylogenetic trees recalled in §2.2. The basis of this 
method is the representation of a network by means of the tripartitions associated 
to its arcs. For each arc e = (u, v) of a DAG N labeled in S, the tripartition of S 
associated to e is 

e{e) = {Al{v), Cl{v) \ Al{v), S \ Cl{v)), 

where moreover each leaf s in Al{v) and Cl{v) \ Al{v) is weighted with the 
greatest number of hybrid nodes contained in a path from v to s (including v and 
s themselves).^ Let 6[N) denote the set of all these tripartitions of arcs of A''. 

In some of the aforementioned papers, the authors enrich these tripartitions 
with an extra piece of information. Namely, they define the reticulation scenario 
RS{v) of a hybrid node v with parents ui,U2 as the set of clusters of its parents: 

RS{v) = {Cl{ui),Cl{u2)}. 

Then, the enriched tripartition ^{e) associated to an arc e is defined as 6{e) if e 
is a tree arc, and as the pair {9{e), RS{v)) if e is a hybridization arc with head v. 
Let W'{N) denote the set of all these enriched tripartitions. 

Actually, Moret, Nakhleh, Warnow et al consider also other variants of this definition, weighting 

only the non-strict descendant leaves or not weighting any leaf, but for the sake of brevity and 
generality we only recall here the most general version. 
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For every T = 9,^, the error, or tripartition, metric relative to Y between two 
DAGs A^i = {Vi,Ei) and N2 = (^2,-^2) labeled in the same set S is defined by 
these authors as 



Unfortunately, and despite the word 'metric', this formula does not satisfy the 
separation axiom on any of the subclasses of phylogenetic networks where it is 
claimed to do so by the authors, and hence it does not define a distance on them: 
for instance, does not satisfy the separation axiom on the class of tree-sibhng 
model phylogenetic networks recalled above. See [7] for a detailed discussion of this 
issue. 

Two other dissimilarity measures considered in [27, 28, 30] are based on the 
representation of a rooted DAG by means of its induced subtrees: the phylogenetic 
trees with the same root and the same leaves as the network that are obtained by 
taking a spanning subtree of the network and then contracting elementary paths 
into nodes. For every rooted DAG A'', let T{N) denote the set of all its induced 
subtrees, and C{N) the set of all clusters of nodes of these induced subtrees. 

Then, for every two rooted DAGs Ni = {Vi,Ei) and N2 = {V2,E2) labeled in 
the same set S, the authors define: 

— m*^^^(A^i, A^2) as the weight of a minimum weight edge cover of the complete 
bipartite graph with nodes T{Ni) U T{N2) and edge weights the value of the 
Robinson-Foulds metric between the pairs of induced subtrees of A^i and ^"2 
connected by each edge. 

— m^'P{Ni, N2) as mr, replacing T by C: 



These measures do not satisfy the separation axiom on the class of tree-sibling 
phylogenetic networks: see, for instance, [27, Fig. 6.8]. On the positive side, Nakhleh 
et al prove in [27, §6.4] and [28, §5] that they are distances on the subclass of time- 
consistent binary galled trees. But it can be easily checked that on arbitrary galled 
trees they do not define distances either: see, for instance. Fig. 3. 

3 Tree-child phylogenetic networks 

Since in this paper we are not interested in the reconstruction of networks, for the 
sake of generality we assume the most general notion of phylogenetic network on 
a set S of taxa: any rooted S'-DAG. So, its hybrid nodes can have any in-degree 
greater than one and any out-degree, and its tree-nodes can have any out-degree. 
In particular, they may contain hybrid leaves and out-degree 1 tree nodes. 

We shall introduce two comparison methods on a specific subclass of such 
networks. 




l/ir(Ari)\r(Ar2)| , \T{N2)\T{Ni)\ 
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Definition 1. A phylogenetic network satisfies the tree-child condition, or it is 
a tree-child phylogenetic network, when every internal node has at least one tree 
child. 

Tree-child phylogenetic networks can be understood thus as general models of 
reticulated evolution where every species other that the extant ones, represented 
by the leaves, has some descendant through mutation. This slightly strengthens the 
condition imposed on phylogenetic networks in [22] , where tree nodes had to have 
at least one tree child, because we also require internal hybrid nodes to have some 
tree child. So, if hybrid nodes are further imposed to have exactly one child (as for 
instance in the definition of model phylogenetic network recalled in §2.3), this node 
must be a tree node: this corresponds to the interpretation of hybrid nodes not as 
individuals but as recombination events, producing a new individual represented 
by their only child. On the other hand, if hybrid nodes represent individuals, then 
a hybrid node with all its children hybrid corresponds to a hybrid individual that 
hybridizes before undergoing a speciation event, a scenario, according to [22], that 
"almost never arises in reality." 

The following result gives two other alternative characterizations of tree-child 
phylogenetic networks in terms of their strict and tree clusters. 

Lemma 2. The following three conditions are equivalent for every phylogenetic 
network N = {V,E): 

(a) N is tree- child. 

(b) Tl{v) / for every node v E V \ Vl. 

(c) Al{v) 7^ for every node v & 

Proof. (a)=^>(b): Given any node v other than a leaf, we can construct a tree path 
by successively taking tree children. This path must necessarily end in a leaf that, 
by definition, belongs to Tl{v). 

(b) ^(c): If f ^ Vl, then, by Corollary 1, 7^ Tl{v) C Al{v), while if u e Vl, 
then, by definition, v G Al{v). 

(c) =>(a): Let v be any internal node. We want to prove that if Al{v) 7^ 0, 
then V has a tree child. So, let s G Al{v)^ and consider the set W of children of v 
that are ancestors of s: it is non-empty, because s must be a descendant of some 
child of V. Let w be a maximal element of W with respect to the path ordering on 
N. If is a tree node, we are done. Otherwise, let v' be a parent of w different 
from V. Let r^v' be any path from a root r to v' . Concatenating this path with 
the arc {v',w) and any path w s, we get a path r s. Since s £ Al(v), this 
path must contain v, and then, since N is acyclic, v must be contained in the path 
r v'. Let w' be the node that follows v in this path. This node w' is a child of 
V and there exists a non-trivial path w' -^w (through v'), which makes w' also an 
ancestor of s. But then w' E W and w' > w, which contradicts the maximality 
assumption on w. □ 
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Next lemma shows that tree-child phylogenetic networks are a more general 
model of evolution under mutation and recombination than the galled trees. 



Lemma 3. Every rooted galled tree is a tree-child phylogenetic network. 

Proof. Let N = (V, E) be a galled tree. If N does not satisfy the tree-child condi- 
tion, then it contains an internal node u with all its children vi,...,Vk G V 
hybrid. 

The node u cannot be hybrid, because in galled trees a hybrid node cannot have 
any hybrid children. Indeed, assume that u has two parents a, 6, and let u' be the 
other parent of the child vi of v. Let x be the least common ancestor of a and h, and 
y the least common ancestor of h and u' . Then the recombination cycles defined by 
the paths (x, . . . , a, u) and (x, . . . , 6, n), on the one hand, and {y, . . . ,b, n, vi) and 
(y, . . . , u' , vi) on the other hand, share the arc (6, u), contradicting the hypothesis 
that is a galled tree. See Fig. l.(a). 

Thus, M is a tree node. In this case, k ^ 2, because galled trees cannot have 
out-degree 1 tree nodes. Now, if u is the root of N, then Al{u) = Vl / and hence, 
by the proof of the implication (c)^(a) in Lemma 2, it has some tree child. If, on 
the contrary, u is not the root of A'', let w be its parent and ui and U2 the parents 
other than u of vi and V2, respectively. Let xi be the least common ancestor of w 
and Ml, and X2 the least common ancestor of w and U2. Then the recombination 
cycles defined by the paths (xi, . . . , ui, f i) and (xi, . . . ,w,u, vi), on the one hand, 
and (x2, . • • , U2,V2) and (x2, . . . ^w^u, V2), on the other hand, share the arc {w, u), 
contradicting again the hypothesis that is a galled tree. See Fig. l.(b). □ 



Fig. 1. (a) In a galled tree, a hybrid node cannot have a hybrid child, (b) In a galled tree, a 
non-root tree node cannot have two hybrid children. 



Remark 1. Not every tree-child phylogenetic network is a galled tree: see, for in- 
stance, the tree-child phylogenetic network in Fig. 4. 



lca{w, ui) 



lca{w, U2) 




10 



We provide now some upper bounds on the number of nodes in a tree-child 
phylogenetic network. 

Proposition 1. Let N = {V, E) be a tree-child phylogenetic network with n leaves, 
(a) \VH\^n-l. 

(h) If N has no out- degree 1 tree node, then |y| ^ 2n — 1 + di{v). 

(c) If N has no out-degree 1 tree node and if m = max{dj(u) | v G Vh}, then 
\V\ ^ (m + 2)(n-l) + l. 

Proof, (a) Let r be the root of N. Consider a mapping t : V \Vl ^ Vt \ {r} 
that assigns to every internal node one of its tree children; since tree nodes have 
a single parent, this mapping is injective. Then, \V\ — \Vl\ ^ \Vt\ — 1 and, since 
= \Vh\ + \Vt\ and \Vl\ = n, 

|y^| = |y|-|VT|^|VL|-l = n-l. 

(b) For every j ^ 2, let Vhj be the set of hybrid nodes with in-dcgree j. If, 
for every hybrid node v, we remove from N a set of di{v) — 1 arcs with head v, 
we obtain a tree with set of nodes V and set of leaves Vj^ (no internal node of 
N becomes a leaf, because when we remove an arc e, since it is a hybridization 
arc, there still remains some tree arc with the same tail as e). Now, in this tree 
there will be at most di{v) nodes with in and out-degree 1: did not have 
any out-degree 1 tree node and, in the worst case, when we remove the di{v) — 1 
arcs with head v, this node and the tails of the removed arcs become nodes of in 
and out-degree 1. Since, in a tree, the number of nodes is smaller than twice the 
number of leaves plus the number of nodes with in and out-degree 1 , the inequality 
in the statement follows. 

(c) If m = m.aji{di{v) \ v G Vh}, then YIv^Vh ^ii'") ^ '^l^ifl- Then, combining 
(a) and (b), 

\V\ ^ 2n-l-h ^ di{v) ^ 2n-l + m\VH\ ^ 2n-l-hm(n-l) = {m+2){n-l)-hl, 
vGVh 

as we claimed. □ 

The upper bounds in Lemma 1 are sharp, as there exist tree-child phylogenetic 
networks for which these inequalities are equalities: for n = 1, point (a) entails 
that is a tree, and (b) and (c) then simply say that N consists only of one node; 
for n ^ 2, see the next example. In particular, for every number n ^ 2 of leaves, 
there exist arbitrarily large tree-child phylogenetic networks without out-degree 1 
tree nodes with n leaves. Of course, if we do not forbid out-degree 1 tree nodes, 
then there exists no upper bound on the number of nodes of the network. 
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Example 1. Let T be the 'comb-like' binary phylogenetic tree labeled in {1, . . . , n} 
described by the Newick string 

(1,(2,(3,. ..,(n-l,n). ..))), 

and let us fix a positive integer number m ^ 2. 

For every i = 1, . . . , n — 1 let us call Vi^m the parent of the leaf i: to simplify 
the language, set Vn,m = n. Notice that vi^m is the root of the tree. Now, for every 
i = l,...,n — 1, split the arc {vi^rn-, i) into a path of length m, 

split the arc (fi,m, "f^i+i,™) into a path of length 2, 

and, for every z = 1, . . . , n — 1 and j = 1, . . . , m — 1, add an arc (wjj, /ii+i). Fig. 2 
displays^ this construction for n = 4 and m = 3. 

The original binary tree had 2n — 1 nodes, and we have added (m — l)(n — 1) 
new tree nodes and n — 1 hybrid nodes (of in-degree m). Therefore, the resulting 
tree-child phylogenetic network has (m -|- 2)(n — 1) -|- 1 nodes. 




Fig. 2. A tree-child phylogenetic network with 4 leaves and 5-3 + 1 nodes. 

In the next sections, we define a distance on the class of all tree-child phy- 
logenetic networks. It is convenient thus to remember here that the tripartition 

^ Henceforth, in graphical representations of phylogenetic networks, and of DAGs in general, 
hybrid nodes are represented by squares and tree nodes by circles. 
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metrics mg or recalled in §2.4 do not define a distance on this class, because 
there exist pairs of non-isomorphic tree-child phylogenetic networks on the same 
set of taxa with the same sets of enriched tripartitions: for instance, the networks 
depicted in Figs. 4 and 8 below (see [7] for details). As far as the metrics m^^'^^ and 
m'^P also recalled in §2.4 goes, they do not define either distances on the class of all 
tree-child phylogenetic networks, because there also exist pairs of non-isomorphic 
tree-child phylogenetic networks on the same set of taxa with the same sets of 
induced subtrees. For instance, the tree and the galled tree depicted in Fig. 3 have 
the same sets of induced subtrees, namely the tree itself, and hence the same sets 
of clusters of induced subtrees. 





Fig. 3. A tree and a galled tree with the same sets of induced subtrees. 



4 The //-representation of tree-child phylogenetic networks 

Let us fix henceforth a set of labels S = {li, . . . , In}: unless otherwise stated, all 
DAGs appearing henceforth are assumed to be labeled in S, usually without any 
further notice. 

Let = (V, E) be an S'-DAG. For every node u £ V and for every i = 1, . . . , n, 
we denote by mi{u) the number of different paths from u to the leaf li. We define 
the path-multiplicity vector, or simply ^-vector for short, of u E F as 

p.{u) = (mi(n), . . . ,m„(u)); 

that is, /x(n) is the n-tuple holding the number of paths from u to each leaf of the 
graph. 

To simplify the notations, we shall denote henceforth by ^^"^ the unit vector 

(0,. .. ,0,1,0,. ..,0). 

^ V ' 

n 

Lemma 4. Let u £ V be any node of an S-DAG N = (V, E). 
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(a) Ifu = lie Vl, then fi{u) = 5f'. 

(b) If u ^ Vl and child(ti) = {vi, . . . ,Vk}, then n{u) = /u(fi) + • • • + //(t'fc). 

Proof. The statement for leaves is triviaL When u G V\Vl, by deleting or prepend- 
ing u we get, for every i = 1, . . . , ra, a bijection 



Remark 2. If u G child(u), then = if and only if v is the only child of u: 
any other child would contribute something to iJi{u). 

Lemma 4 implies the simple Algorithm 1 to compute the jU- vectors of the nodes 
of an S'-DAG in polynomial time. Since the height of the nodes can be computed in 
0{n+\E\) time, it takes 0{n\E\) time to compute ix{N) on an S'-DAG = {V, E) 
with n leaves. 

Algorithm 1. Given an S-DAG N = {V,E), compute n{N). 

begin 

for i = 1 , . . . , n do 

set = 
sort V \Vl increasingly on height 
for each x e V \ Vl do 

let ui, . . . ,yk £ V be the children of x 
set i^{x) = n{yi) H h n{yk) 



Example 2. Consider the tree-child phylogenetic network depicted in Fig. 4. Ta- 
ble 1 gives the //-vectors of its nodes, sorted increasingly by their heights. 



node 


height 


^-vector 


node 


height 


/i-vector 


node 


height 


/i-vector 


1 





(1,0,0,0,0) 


C 


1 


(0,0,0,1,0) 


a 


6 


(1,1,1,1,0) 


2 





(0,1,0,0,0) 


f 


2 


(0,0,1,1,0) 


d 


6 


(0,1,1,1,1) 


3 





(0,0,1,0,0) 


B 


3 


(0,0,1,1,0) 


c 


7 


(0,1,1,2,1) 


4 





(0,0,0,1,0) 


e 


4 


(0,1,1,1,0) 


b 


8 


(0,1,2,3, 1) 


5 





(0,0,0,0,1) 


A 


5 


(0,1,1,1,0) 


r 


9 


(1,2,3,4, 1) 



Table 1. /j- vectors of the nodes of the network depicted in Fig. 4. 



Exam,ple 3. Consider the phylogenetic network depicted in Fig. 5. Table 2 gives 
the /Lt-vectors of its nodes sorted increasingly by their heights. 




which clearly implies the statement in this case. 



□ 



end 
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Fig. 4. The tree-child phylogenetic network used in Example 2. 



node 


height 


/i-vector 


node 


height 


/^-vector 


node 


height 


/i-vector 


1 





(1,0,0,0,0) 


B 


1 


(0,1,0,0,0) 


/ 


3 


(0,1,1,1,0) 


2 





(0,1,0,0,0) 


C 


1 


(0,0,1,0,0) 


c 


4 


(1,1,1,1,0) 


3 





(0,0,1,0,0) 


D 


1 


(0,0,0,1,0) 


e 


4 


(1,1,1,1,0) 


4 





(0,0,0,1,0) 


a 


2 


(1,1,0,0,0) 


d 


5 


(1,1,1,1,1) 


5 





(0,0,0,0,1) 


9 


2 


(0,0,1,1,0) 


r 


6 


(2,2,2,2,1) 


A 


1 


(1,0,0,0,0) 


b 


3 


(1,1,1,0,0) 





Table 2. /^-vectors of the nodes of the network depicted in Fig. 5. 



Definition 2. The /u-representation of a DAG N = {V,E) is the multiset ^{N) 
of ^-vectors of its nodes: its elements are the vectors fJ-{u) with u £ V, and each 
one appears in /i(iV) as many times as the number of nodes having it as fi-vector. 

It turns out that a tree-child phylogenetic network can be singled out up 
to isomorphism among all tree-child phylogenetic network by means of its fi- 
representation (Thm. 1). Before proceeding with the proof of this fact, we establish 
several auxiliary results. 

The following lemma shows that the path ordering on a tree-child phylogenetic 
network is almost determined by its /^-representation. In it, and henceforth, the 
order ^ considered between /x- vectors is the product partial order on N"". 

{nil, • • • 1 fnn) ^ ("1-1, • • • 1 fn'n) *^==^ ^ "^i for every i = 1, . . . , n. 

Lemma 5. Let N = (y,E) be a tree- child phylogenetic network. 

(a) If there exists a path u-^v, then ^{u) ^ /i(t'). 

(b) If fi{u) > fi{v), then there exists a path u-^v. 

(c) If fi{u) = then u and v are connected by an elementary path. 
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Fig. 5. The tree-sibling, non-tree-child, phylogenetic network used in Example 3. 

Proof. Assertion (a) is a straightforward consequence of Lemma 4. 

As far as (b) and (c) goes, let us assume for the moment that fi{u) ^ n{v) 
and let k G Ti{v)] in particular, mi{y) ^ 1. Then, mi{u) ^ rni{v), and therefore 
there exists also a path u-^li. Now, consider the tree path v-^li. By Lemma 1, 
it must happen that either the path u-^li contains the path v-^li or vice versa, 
and therefore u and v are connected by a path. 

If > by (a) there cannot exist any path v-^u, and therefore there 

exists a path u-^v: this proves (b). 

On the other hand, if fi{u) = fi{v), nothing prevents the existence of a path 
u -w ti or a path v u. To fix ideas, assume that there exists a path u-^ v, say 
(■u, . . .,Vk-i,v). Since, by (a), 

p,{u) ^ ^{vi) ^ • • • ^ ^(^^fc-i) ^ /"(-y) = Ai('u), 

we conclude that 

^{u) = fi{vi) =■■■ = n{vk-i) = 

As we noticed in Remark 2, this implies that each one of n, . . . , v^-i has only one 
child, the node that follows it in this path. Therefore the path u-^v is elementary. 
The same argument shows that if it is the path v u which exists, then it is 
elementary. □ 

Remark 3. Assertion (a) in the previous lemma holds for every DAG, but asser- 
tions (b) and (c) need not be true if the DAG N does not satisfy the tree-child 
condition, even if is a binary tree-sibling phylogenetic network. Indeed, consider 
the phylogenetic network described in Fig. 5. We have in it that 

At(c) = (1,1,1,1,0) =/i(e) 
Md) = (1,1, 1,1,1) 
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Then, c and d have the same /U- vectors but they are not connected by any path, 
yielding a counterexample to (c). And //(d) > //(e), but there is no path d-^ e, 
which yields a counterexample to (b) . 

The next two lemmas show how to recover the children of a node in a tree-child 
phylogenetic network from the knowledge of the jU-representation of the network. 

Lemma 6. Let N be any DAG. Let u G V be any internal node, and let 

Mu = {w eV \ u> w}. 

Then, M„ has maximal elements, and all of them are children of u. 

Proof. The set is non-empty, since u is not a leaf and every descendant of u is 
in M„. Since is finite, it has maximal elements. Let v be any such a maximal 
element. Since u > v, there exists a non trivial path u ^ v. If this path passes 
through some other node w, then u > w > v, against the assumption that v is 
maximal in M^. Therefore, the path u-^v has length 1 and is a child of u. □ 

The maximal elements of are exactly the children of u such that the arc 
{u, v) is the only path u-^ v. This includes all tree children of u, as well as all 
hybrid children v oi u such that no other parent of v is a descendant of u. In time 
consistent phylogenetic networks, this covers all children of u. But in arbitrary 
tree-child phylogenetic networks, this need not cover all of them. Consider for 
instance the right-hand side tree-child phylogenetic network in Figure 3. In it, the 
only maximal element of Mb is a, but A is also a child of 6. 

Lemma 7. Let N be a tree- child phylogenetic network. Let u & V be any internal 
node and vi, . . . ,Vk some of its children. 

(a) If ii{u) = iJ,{vi) -|- • • • -|- ^{vk), then u has no other children. 

(b) If ti{u) > ii{vi) -\ h n{vk), let 

Mu,vi,...,vk = {w eV \ u> w, ii{u) ^ ii{w) + ii{vi) -\ F- Ai(-Ufc)}. 

Then, Mu^vi,...,Vk has maximal elements, and all of them are children of u and 
different from vi, . . . ,Vk. 

Proof. Let us assume that = n^vi) -|- • • • -|- //(ufc). Then u cannot have other 
children, since if it has any other child w, then 

//(n) ^ fi{vi) -\ h fj,{vk) + tJ,{w) > fj,{vi) -\ h fJ-{vk). 

Assume now that iJ,{u) > fi{vi) -I- • • • -I- lJ.{vk). Then, by Lemma 4.(b), u has 
other children than vi,. . . ,Vk. Let N' be the DAG obtained by removing from 
A'' the arcs (n, ui), . . . , (u,Vk). For any node s € V, let m'-(s) be the number of 
different paths s li in A', and set //'(s) = {m[{s) , . . . , m'^{s)) . Then, /Li'(tt) = 
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n{u) — (fJ'ivi) + • • • + fi{vk)), because the paths u U in that are not in N' 
are exactly those whose first visited vertex is one oi vi, . . . ,Vk- Moreover, if w is 
a descendant of u in iV, then ^x'{w) = iJ,{w), because no path w li in N can 
possibly contain any arc {u, Vi) (it would form a cycle with the path u-^w). 
Then, we have that = for every w G Mu,vi,...,vkJ ^'^d thus 

Mu,vi,...,vk = {'wEV\u>whiN and //'(n) ^ n'{w)}. 

Now, it turns out that w £ Mu^vi,...,vk if) only if, there exists a non-trivial 
path u ^ w in A^'. Indeed, if there exists a non-trivial path u w in N' , then 
there exists also the same path in N, and hence u > w in N, and moreover, by 
Lemma 5. (a) above, /u'(n) ^ fi'(w). Conversely, let w be a descendant of u in N 
and assume that fj,'{u) ^ fj'(w) = fi(w). If li G Tl(w), then ra',-{u) ^ mi{w) = 1. 
Take the tree path w^li in A^, which also exists in N' , and any path u-^li in A^'. 
By Lemma 1 (applied to N'), either u-^k contains w-^k or vice versa. But the 
existence of a non-trivial path u-^w in N prevents the existence of a path w-^u 
in N' . Therefore, it is u-^li that contains w-^li and in particular there exists a 
path u-^w also in N' . 
So, 

Mu,vi,...,vk = {weV\u>w in N'}. 

Since at least one child of u has survived in N', the previous lemma implies that 
this set has maximal elements and they arc children of u in N', and hence they 
are also children of u in A'^ and they are different from vi, . . . ,Vk. □ 

As we have already mentioned. Lemma 6 applies to any DAG (and we make 
use of this fact in the proof of Lemma 7, because the DAG N' we consider in 
it need not satisfy the tree-child condition, and can have more than one root as 
well as out-degree 1 tree nodes), but Lemma 7.(b) need not be true if A^ does 
not satisfy the tree-child condition. Consider again, for instance, the tree-sibling 
phylogenetic network A^ described in Figure 5. In it, c is a maximal element of 
Mr.d = {x £ V \ r > X, fi{r) ^ fj,{x) + fi{d)}, but it is not a child of r. 

We can prove now our main result. 

Theorem 1. Let N,N' be tree-child phylogenetic networks. Then, N = N' if, and 
onlyif,ij{N)=fi{N'). 

Proof. Let N = (V, E) be a tree-child phylogenetic network labeled in S, and let 

n{N) be its ^u-representation. Let C x N be the set consisting of the vectors 
of the form (x, i) with x G A*(A/") and i between 1 and the multiplicity of x in fJ,{N). 
Consider on the partial order y defined by 

{x, i) >- {y,j) ■^=^ x > y with respect to the product partial order, 
or: X = y and i < j. 
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Wc know from Lemma 5 that if x G N" belongs to /i(A^) with multiphcity 
m ^ 1, then there exist m nodes in A'^ with //-vector x, and that they form an 
elementary path. For every node v eV, let be the position of v in the elementary 
path formed by all nodes with the same /t-vector as v. In particular, if //(u) appears 
in l-i{N) with multiplicity 1, then = 1 

Lemma 5 implies then that the mapping 

V ^ fl{v) = {fi{v),iy) 

is an isomorphism of partially ordered sets between V with the path ordering and 
with the partial order y. Indeed, ifu > v then either n{u) > ii{y) or = ii{y) 

and u appears before v in the elementary path of all nodes with this /i-vector, and 

hence iu < iy Conversely, if > fi{v) or if //(u) = fi{v) and iu < iv, then there 

exists a non-trivial path u-^v. 

Therefore, we can rephrase the last two lemmas as follows: 

(1) For every u & V internal, the set 

Mu = {fi{w) e fi{y) I fiiu) >- 

has maximal elements, and all of them are images under p, of children of u. 

(2) For every u G V and for every vi, . . . ,Vk G child('u), 

(a) If iJ,{u) > fJ,{vi) + ■ ■ ■ + v{vk), then the set 

Mu,vi,...,vk={P'{'(") e fi{V) I fi{u) >~ fi{w), n{u) ^ ii{w)+ix{vi)-\ \-n{vk)} 

has maximal elements, and all of them are images under Ji of children of 
u other than vi, . . . , Vk- 

(b) If fj,{u) = iJ,{vi) -|- • • • -|- /u('Ufe), then u has no other children. 

Wc shall prove that wc can recover the set E of arcs in N from n{N). To do that, 
consider the set C x V"^ obtained through the application of Algorithm 2: 

Algorithm 2. Given IJ.{N), compute E^. 

begin 

set Ef, = $ 

sort Vn decreasingly on the partial order >z 
for each {x,i) G do 
set m = X 
while m > do 

for each {y,j) G such that {x,i) >- {y,j) do 
if m ^ y then 

add the arc {{x,i), {y,j)) to E^ 
set m = m — y 

end 
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Let us prove that {u,v) G E if, and only if, {fl{u) , fl{v)) G E^. To do that, 
let M G y be an arbitrary node of A''. If there is no arc in E^ with source fi{u), 
it can only be because there is no v G F such that iJ,{u) > n{v), and hence u 
is a leaf of N and E does not contain any arc with source u, either. Otherwise, 
let {jl{u) , Jl{vi)) , {fj,{u) , p,{v2)) , . . . , {fl{u), fi{vi)) be the arcs contained in E^ with 
source fi{u), given in the order they are added to E^. This entails that jl{vi) is a 
maximal element of 

{(l{v) G fl{V) I fi{u) y jl{v)} = Mu 
and, for each i = 2, . . . , Z, lx{vi) is a maximal element of 

{ll{v) G il{V) I Ji{u) y fi{v),n{u) - {n{vi) + ... + l^{vi-i)) ^ i^{v)} _ 

= {v G V I fi{u) y il{v), iJ.{u) ^ iJ,{v) + i2{v-i) + . . . + |J.{vi--l)} = Mu,vl,...,vi-l■ 
ThcIc^o'cc, as we have recalled in points (a) and (b.l) above, the nodes vi, . . . ,vi 
are children of u in N, that is, {u, vi), {u,V2), ■ ■ ■ , {u,vi) G E. On the other hand, 
the algorithm adds arcs {li{u) , p,{vi)) to E^^ until it happens either that fi{u) = 
+ • • • + f^{vi), in which case, by (b.2), vi, . . . ,vi are exactly the children of 
u in N, or that the set of nodes is exhausted and //(u) > n{vi) + • • • + //(w;): 
but the latter cannot happen, because fi{u) must be the sum of the /x-vectors 
of its children in A^. Thus, in summary, child(«) = {vi, . . . ,vi} in N and thus 
{u, vi), {u, V2), . . . ,{u, vi) are also all the arcs contained in E with source u. 

This proves that E = {{u,v) | p,{u),jl{v)) G -E^}, as we claimed. Now, if 
N = {V,E) and N' = {V',E') are two tree-child phylogenetic networks such that 
fi{N) = n{N'), we have bijections 

V ^v^ = v;,^ v. 

Renaming in V the nodes of N' through this bijection V — > V, we obtain a 
phylogenetic network A^" = (V, E") isomorphic to A^' and such that ijl{N") = 
IJl{N') = fi{N). Let Efj, and E'^ be, respectively, the sets or arcs obtained by 
applying the previous algorithm to N and N". Since n{N) = iJ,{N"), we have that 
Eij, = E'^, and hence 

E = {{u,v) I fliu),fliv)) G E^} = {{u,v) I fliu),fliv)) G 4'} = E". 

This implies that A^ = A^" and therefore N = N' as DAGs. Now, this isomorphism 
clearly preserves the leaves' labels, because it preserves /x-vectors. Therefore, N = 
N' also as S-DAGs. 

This proves the "if" implication in the statement. Of course, the "only if" 
implication is obvious. □ 

To recover, up to isomorphism, a tree-child phylogenetic network A^ from its /i- 
representation n{N), it is enough to compute the set associated to the multiset 
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l-i^N), then to apply Algorithm 2 to compute the set of arcs E'^, and finally to 
label each leaf of the resulting DAG, which will have the form ((5|"\mj) with rrii 
the multiplicity of J^"^ in n{N), with the corresponding label li. 

Example 4- Let us apply this procedure to the /^-representation of the tree-child 
phylogenetic network N depicted in Fig. 4. From the multiset IJ,{N) described in 
Table 1, we obtain the following set V^, which we give sorted decreasingly on >- 
(and, to simplify the description of the application of the algorithm, we give names 
Xi to its elements): 

V;, = {xi = ((1,2,3,4,1),!), X2 = ((0,l,2,3,l),l), = ((0, 1, 1, 2, 1), 1), 

X4 = ((0,ia,l,l),l), a;5 = ((1,1,1,1,0),!), xg = ((0, 1, 1, 1, 0), 1), 
X7 = ((0,1, 1,1,0), 2), X8 = ((0,0,1, 1,0),!), X9 = ((0,0,1, 1,0), 2), 
xio = ((1,0,0,0,0),!), xii = ((0,1,0,0,0),!), X12 = ((0, 0, 1, 0, 0), 1), 
xi3 = ((0,0,0,1,0),!), xi4 = ((0,0,0,!,0),2), X15 = ((0,0,0,0, 1), !)}. 

We shall denote the first and the second component of each Xk by and ik, 
respectively. 

We begin with an empty set of arcs: 



Then we proceed with the for each in Algorithm 2, visiting all elements of in 
the given order: 

(xi) We set m = (1,2,3,4,1). Then, since m ^ iJ,2, we add (xi,X2) to Efj_ and we 
set m = m — 1^2 = (1)1)1)1)0). The first next element Xk of with /i/fc ^ m is 
X5. Then, we add (xijXs) to and we set m = m — /X5 = (0,0,0,0,0). This 
makes us to stop with xi. At the end of this step, we have 

= {(xi,X2),(xi,X5)}. 

(X2) We set m = (0,1,2,3,1). Since m ^ ^^3, we add (x2,X3) to i?^ and we set 
m = m — Us = {0,0, 1, 1, 0). The first next element Xfc of with /Hfc ^ m is xg. 
Then, we add (x2, xs) to E^, we set m = m — fig = (0, 0, 0, 0, 0), and we stop. 
So, at the end of this step, we have 

E^t = {(X1,X2),(X1,X5),(X2,X3),(X2,X8)}, 

(xs) We set m = (0,1,1,2,1). Then, since m ^ ;U4, we add (x3,X4) to and we 
set m = m — /X4 = (0, 0, 0, 1, 0). The first next element Xk of with ^ m is 
X13. Then, we add (x3,xi3) to we set m = m — /X13 = (0,0,0,0,0), and we 
stop. At the end of this step, we have 

E^, = {(Xi, X2), (Xi, X5), (X2, X3), (X2, Xs), (X3, X4), (X3, X13)}. 
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(x4) We set m = (0, 1, 1, 1, 1). The first element Xk with ^ m is xq, and therefore 
we add (x4, xq) to and we set m = m — //e = (0, 0, 0, 0, 1). The first next 
element x^ of with /Ufe ^ m is X15. Then, we add (x^jXi^) to E'^u, we set 
m = m — )Ui5 = (0, 0, 0, 0, 0), and we stop. At the end of this step, we have 

E^, = {{xi,X2), ixi,X5), (X2, X3), (X2, Xg), (xs, X4), {X3, X13), (X4, Xe), (X4, Xis)}. 

(X5) We set m = (1,1,1,1,0). Since m ^ ^ug, we add (xsjXg) to i?^ and wc set 
m = m — jiQ = (1,0,0,0,0). The first next element x^ of with /i/fe ^ m is 
xiQ. Then, we add (x5,a:io) to E^^^. Now m = m — xio = (0,0,0,0,0), and we 
stop. At the end of this step, we have 

Ef, = {{X1,X2), {X1,X5), {X2,X3), {X2,X8), {X3,X4:), (a;3,Xl3), (X4,X6), (a;4,Xl5), 
{x5,Xq),{x5,Xio)}. 

(xq) We set m = (0,1,1,1,0). Since m ^ /ny, we add {xq,X7) to E"^ and we set 
m = m — fi7 = {0,0, 0, 0, 0), and we stop. At the end of this step, we have 

E^, = {{XI,X2), {XI,X5), {X2,X3), {X2,X8), (a;3,X4), (X3,X13), {x4,Xg), (x4,X15), 
{x5,Xq),{x5,Xio),{X6,X7)}. 

(xy) We set (again) m = (0, 1, 1, 1, 0). Since m ^ ns, we add {xj, xg) to E^ and we 
set m = m — /X8 = (0, 1,0,0,0). The first next element Xk of with /Xfe ^ m is 
xu- Then, we add (,T7,,tii) to i?^. Now m = m — /xn = (0,0,0,0,0), and we 
stop. At the end of this step, we have 

En = {{xi,X2), (Xi, X5), (X2, X3), (X2, Xs), {xs, X4), (^3, X13), {X4, Xq), (x4, X15), 
(X5, Xq), (x5, Xio), (x6, X7), (X7, Xg), (X7, Xii)}. 

(xs) We set m = (0,0,1,1,0). Since m ^ fig, we add {x8,xg) to ^ji^, we set m = 
m — jUg = (0, 0, 0, 0, 0), and we stop. At the end of this step, we have 

En = {(Xi, X2), (xi, X5), (X2, X3), (X2, Xg), (X3, X4), (X3, X13), (x4, Xe), (X4, X15), 
(X5, Xe), (X5, Xio), (X6, X7), (X7, xs), (X7, Xll), (X8, Xg)}. 

(xg) We set m = (0, 0, 1, 1, 0). The first element Xk in with < m is X12, and 
then we add (xg, X12) to E^ and we set m = m — /X12 = (0, 0, 0, 1,0). The first 
next element x^. in with ^ m is X13. Then, we add (x9,xi3) to E^,- Now 
m = m — /X13 = (0, 0, 0, 0, 0), and we stop. At the end of this step, we have 

En={{xi,X2), (Xi, X5), (X2, X3), (X2, Xs), (X3, X4), (X3, X13), (X4, Xq), (x4, X15), 

(X5, Xe), (X5, Xio), (X6, X7), (X7, Xs), (X7, Xll), (xs, Xg), (xg, X12), (xg, X13)}. 

(xio) Since m = (1, 0, 0, 0, 0) has only one non-zero entry, and xio is the only element 
of with /ifc = m, the algorithm does not find any arc in E^ with head xiq. 
The same happens with xn, xi2, and X15. 
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(xis) We set m = (0,0,0,1,0). Since m ^ /ii4, we add (xi3,xi4) to we set 
m = m — /ii4 = (0, 0, 0, 0, 0), and we stop. At the end of this step, we have 

E^_i={{xi,X2), (a;i,X5), (X2,X3), {X2,X^), (X3,X4), (xSjXis), (X4,X6), (x4,Xi5), 

(x5, xe), (x5, xio), (a::6, a;?), (a;?, a^s), [xi, xn), (xs, xg), (a^g, a;i2), (a;g, xiz), 
(a;i3,a;i4)}- 

(X14) Since m = (0,0,0, 1,0) has only one non-zero entry, and no other element of 
Vfj_ after a;i4 has this first component, the algorithm does not find any arc in 
with head xn. 

The DAG (V^, E^) obtained up to now is depicted in Fig. 6. Finally, we would 
label 1, 2, 3, 4, and 5 the nodes xio, xn, X12, X14, and xis, respectively. The result- 
ing DAG labeled in {1, ... ,5} is clearly isomorphic to the tree-child phylogenetic 
network in Fig. 4. 




Fig. 6. The DAG recovered from the ^-representation of the phylogenetic network in Fig. 4. 

Remark 4- The thesis of Theorem 1 need not hold if A'^i and A^2 do not satisfy 
the tree-child condition. Indeed, it is not difficult to check that the tree-sibling 
phylogenetic network given in Figure 7 has the same /i-representation as the one 
given in Figure 5, but they are not isomorphic as S'-DAGs. 

5 The /^-distance for tree-child phylogenetic networks 

For every pair of DAGs Ni and N2 labeled in the same set S, let 
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Fig. 7. This phylogenetic network has the same /^-representation as the one in Fig. 5. 

where the symmetric difference A refers to multisets: if a vector belongs to /u(A'^i) 
with multiplicity a and to f^{N2) with multiplicity b, then it belongs to ^(A^i) A 
fj,{N2) with multiplicity |a — 6|, and hence it contributes |a — 6| to \fi{Ni) A/i(iV2)|. 

Theorem 2. Let Ni, N2, be tree-child phylogenetic networks on the same set 
of taxa. Then: 

(a) d^{Ni,N2)^0 

(b) d^{Ni,N2) = if, and only if Ni ^ N2 

(c) d^{Ni,N2) = d^{N2,N,) 

(d) d^{Ni,N3) ^ d^{Ni,N2)+d^{N2,Ns). 

Proof, (a), (c), and (d) are direct consequences of the properties of the symmetric 
difference, and (b) is a consequence of Theorem 1. □ 

Therefore, defines a distance on the class of all tree-child phylogenetic net- 
works: we shall call it the ^-distance. 

We have shown in Sect. 4 that the //-representation of an S-DAG can be com- 
puted in polynomial time. Now, given two 5-DAGs A^i = {Vi,Ei) and N2 = 
(^2,-^2) with n leaves together with their //-representations /t(iVi) and /i(iV2), the 
simple Algorithm 3 performs a simultaneous traversal of the internal nodes of A'^i 
and A^2j sorted by their /i-vectors, in order to compute the /t-distance d^{Ni,N2) 
in 0(n|y|) time, where \ V\ = max(|yi|, IV2I). 

Algorithm 3. Given the ^-representations n{Ni) and fi{N2) of two S-DAGs Ni = 
{Vi.Ei) and N2 = {V2,E2), compute d^{Ni,N2). 

begin 

let Vl be the set of leaves, common to A'^i and A'^2 
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sort Xi = Vi \ iVijL increasingly according to the lex ordering of the ^u- vectors 
sort X2 = V2 \ (V2)l increasingly according to the lex ordering of the //-vectors 
set d = 

while Xi / and X2 / do 

let xi and X2 be the first element of X\ and X2, respectively 
case ii{xi) < ii{x2) 

set Xi = Xi\ {xi} 

set d = d + 1 
case ii{xi) > /i(x2) 

set X2 = X2 \ {X2} 

set d = d + 1 
otherwise 

set Xi = Xi\ {xi} 

set X2 = X2 \ {x2} 
return d+ \Xi\ + IX2I 
end 

Example 5. Consider the tree T and the network N in Fig. 3. Their /x-representa- 
tions are 

M(r) = {(1,0,0), (0,1,0), (0,0,1), (0,1,1), (1,1,1)} 

li{N) = {(1, 0, 0), (0, 1, 0), (0, 0, 1), (0, 1,0), (0, 1, 1), (0, 2, 1), (1, 2, 1)} 

and therefore 

/x(T) A fi{N) = {(0, 1,0), (1, 1, 1), (0, 2, 1), (1,2, 1)} 
from which we obtain d^{T, N) = 4. 

Exam,ple 6. It was shown in [7] that the tree-child phylogenetic network N depicted 
in Fig. 4 could not be distinguished from the tree-child phylogenetic network N' 
depicted in Fig. 8 below using the tripartition metric. We have already given in 
Table 1 the /U-vectors of the nodes of N . In Table 3 we give the /x-vectors of the 
nodes of N' . From these tables we get that 

iM(iV) A /x(7V')| = {(0, 1, 1,2, 1), (0, 1,2, 2, 1)}, 

which implies that d^^N, N') = 2. 

Example 7. Let N be any tree-child phylogenetic network labeled in S and let v 
be any internal node of it. Let N' be a tree-child phylogenetic network obtained 
by adding to A'^ a new internal node v', an arc {v,v'), and then distributing the 
children of v between v and v' so that N' remains tree-child and v' does not become 
a leaf. Then //(iV) A fi{N') = {n{v')} and therefore d^{N, N') = I. 
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Fig. 8. The tree-child phylogenetic network A'^' compared in Example 6 with the tree-child phy- 
logenetic network A'^ from Fig. 4. 
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Table 3. p-vectors of the nodes of the network depicted in Fig. 8. 



So, expanding a node into an arc yields //-distance 1, just as it happens with 
Robinson- Foulds distance for phylogenetic trees. This is consistent with the fact, 
which we shall prove later, that the //-distance extends the Robinson-Foulds dis- 
tance to tree-child networks: cf. Theorem 3 below. But, contrary to the tree case, 
two tree-child phylogenetic networks can be at //-distance 1 without any one of 
them being obtained by expanding a node into an arc in the other one. Consider, 
for instance, the tree-child phylogenetic networks and A^' labeled in {1,2,3} 
depicted in Fig. 9. Their //-representations are 

fi{N) = {(1,0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 0), (1, 1, 0), (1, 1, 1), (2,2, 1)} 
fi{N') = {(1,0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 0), (1, 1, 0), (2, 2, 1)} 

and thus df,{N,N') = 1. 

It should also be noticed that, also against what happens in the tree case, 
collapsing an arc in a tree-child phylogenetic network N into a node (that is, given 
an arc {v,v'), removing v' and this arc, and replacing every other arc with tail 
or head v' by a new arc with tail or head, respectively, v) need not produce a 
network at /i-distance 1 of A'^: for instance, if v and v' hybridize in N, or if v' is 
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a non- strict descendant of v. We leave to the interested reader to draw specific 
counterexamples . 



Fig. 9. Two tree-child phylogenetic networks A'' (left) and N' (right) at /x-distance 1. 



Example 8. There exist 66 pairwise non-isomorphic binary tree-child phylogenetic 
networks with 3 leaves. All of them have an even number of internal nodes, and 
therefore the ;U-distance between two of them is always an even number. In Propo- 
sition 2 below we shall see that this /i-distance is smaller than or equal to 12. The 
left-hand side histogram in Fig. 10 shows the distribution of distances between 
unordered pairs of such networks. 

In a similar way, there exist 4059 pairwise non-isomorphic binary tree-child 
phylogenetic networks with 4 leaves. Again, all of them have an even number of 
internal nodes, and therefore the //-distance between two of them is always even, 
and in Proposition 2 we shall see that it is smaller than or equal to 18. The right- 
hand side histogram in Fig. 10 shows the corresponding distribution of distances. 

See the Supplementary Material for more details. 

Every phylogenetic tree is a tree-child phylogenetic network, and, as we have 
already mentioned, it turns out that the restriction of this //-distance to the class 
of phylogenetic trees is the Robinson-Foulds metric. 

Theorem 3. For every phylogenetic trees Ti,T2 on the same set of taxa S, 



Proof. Let S = {li, . . . , In}- The uniqueness of paths in trees implies that, if T = 
{V, E) is any phylogenetic tree labeled in 5, then, for every tt E F, 




dt,{Ti,T2) = dRF{Ti,T2). 
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Fig. 10. Histograms of /i-distances between unordered pairs of binary tree-child phylogenetic trees 
with 3 (left) and 4 (right) leaves. 



Therefore, the //-vector of a node « of a phylogenetic tree labeled in S is the image 
of its cluster Cl{u) under the bijection between the powerset V{S) of S and {0, 1}" 
that sends each subset ^ of 5 to its characteristic vector 

XA = iXAih), ■ ■ .,XA{ln)), with XA{k) = 

Then, given two phylogenetic trees Ti and T2 on the set of taxa S, this bijection 
V{S) — > {0, 1}" transforms the sets Cl(Ti) and Cl(T2) of clusters of their nodes 
into their //-representations /t(Ti) and /((T2), respectively, and hence the symmetric 
difference of the former into the symmetric difference of the latter. Therefore, 

|7r(ri) A 7r(r2)| = \Cl{T^) A Cl{T2)\ = |//(Ti) A //(Ta)!, 

as we claimed. □ 

The //-distance takes integer values. Its smallest non-zero value is 1 but it 
can be arbitrarily large. If we bound the in-degree of the hybrid nodes of the net- 
works, then we can compute the diameter of the resulting subclass of phylogenetic 
networks. 

Let TCAA"'™" be henceforth the class of all tree-child phylogenetic networks on 
a fixed set of taxa of n elements, without out-degree 1 tree nodes and with all their 
hybrid nodes of in-degree at most m. 

Proposition 2. For every Ni,N2€ TCM'^''^, 

df,{Ni,N2) ^ 2(m + l)(n-l) 

and there exist pairs of networks in TCN^'^ at ^-distance 2{m + l)(n — 1). 

Proof. Let Ni, N2 G TCA^"'™. By Proposition l.(c), each one of /t(iVi), /i(Af2) 
has at most (m -|- 2)(n — 1) -|- 1 elements, from which at least the n /i- vectors 
corresponding to the leaves will appear in both sets. Therefore 

di_,{Ni,N2) ^2((m + 2)(n-l) + l) - 2n = 2(m + l)(n - 1). 



r 1 like A 
I0 iili^A 
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To find a pair of networks in TCAA"'™ at distance 2(m + 1) (n — 1), let be the 
tree-child phylogenetic network with n leaves and n — 1 hybrid nodes of in-degree 
m described in Example 1. A simple argument by induction shows that 



H{hn) 


= (0,. 


.,0,0,0,0,1) 




= (0,. 


. , 0, 0, 0, 1, k) for every k = 1, . . . ,m 


/U(/in-l) 


= (0,. 


., 0,0,0,1, m) 




= (0,. 


. , 0, 0, 1, fc, km) for every k = 1, . . . ,m 


IJi{hn-2) 


= (0,. 


. ,0, 0, l,m, m^) 


KVn-3,k) 


= (0,. 


. , 0, 1, k, km, km^) for every k = l,... 



and, in general, 

n-j-l 

m^) for every j = 0, . . . , n — 2 
km^) for every k = 1, . . . ,m, j = 1, . . . ,n— 1 

Therefore, fJ-{N) contains, beside the /^-vectors 6^^^ of leaves, all vectors of the form 

n-j-2 

(0, . . . , 0, 1, k, km, . . . , km^), k = 1, . . . ,m — 1, j = 1, . . . ,n — 1 
with multiplicity 1 , and all vectors of the form 

n-j-l 

{0,...,0,l,m,m'^,...,m^), j = l,...,n-2 
with multiplicity 2. 

Now let A^' be the tree-child phylogenetic network in TCM"'''^ obtained by per- 
forming the same construction starting with the binary phylogenetic tree described 
by the Newick string 

in, (n-l,(n-2,...,(2,l)...))). 

The same argument shows that /x(A^') contains, again beside the /x-vectors 6^^^ of 
leaves, all vectors of the form 

n-j-2 

{km-' , . . . , km, A;, 1, 0, . . . , 0), A; = 1, . . . , m — 1, j = 1, . . . , n — 1 

with multiplicity 1 , and all vectors of the form 

n-j-l 

{m^,...,m^,m,l,Q,...,Qi), j = l,...,n-2 



li{hn-j) = (b, . . . , 0, 1, m, m^, . . . , 

n-j-2 

IJ-{vn-j,k) = (0, . . . , 0, 1, fc, km, 
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with multiplicity 2. 

Then, IJ,{N) and n{N') have no //-vector of internal node in common, and 
since each one has (m + l)(n — 1) internal nodes, this implies that d^(iV, iV') = 
2(m+l)(n-l). □ 

This result allows us to normalize the /x-distance on TCA/""'"*. 

Corollary 2. The mapping 

" 2(^+l)(n-l) '''''^'-^^' 

is a distance TCAA"'™ i/iai taA;es values in the unit interval [0, 1]. 



6 The alignment of tree-child phylogenetic networks 

Let A^i = {Vi,Ei) and A''2 = (V2,-£'2) be two tree-child phylogenetic networks 
labeled in the same set 5 = {Zi, . . . , In}- For simplicity, we assume that they don't 
have out-degree 1 tree nodes, and therefore, if two nodes in one of these networks 
have the same /x-vector, then they must be a hybrid node and its only child. 
For every vi € V\ and V2 let 

n 

\mi{vi) - mi{v2)\ 

i=l 

Jo if vi,V2 are of the same type (both tree nodes or both hybrid) 
\ 1 if vi,V2 are of different type 

Notice that H{vi,V2) is the Manhattan, or Li, distance between ^i{t)i) and n{v2). 
The advantage of this distance over the Euclidean distance is that it takes integer 
values on N'*. 

Define finally the weight of the pair {vi,V2) as 

W{Vl,V2) = H{Vl,V2) H — — . 

To fix ideas, assume that |Fi| ^ \V2\. Then, given a matching between Ni and 
N2, that is, an injective mapping M : V\ ^ V2 that preserves leaves and their 
labels, its total weight is defined as 

w{M) = ^ w{v,M{v)). 
veVi 

An optim,al alignment between A'^i and is a matching with the smallest total 
weight. Such an optimal alignment can be computed in time 0((|Vi|-|-|l^|)^) using 
the Hungarian algorithm [17,24]. 



H{vi,V2) = 
X{vi,V2) = 
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Proposition 3. A matching M between Ni and N2 is an optimal alignment if, 
and only if, it minimizes the sum 

J2 H{v,M{v)) 
vGViXVl 

and, among those matchings minimizing this sum, it maximizes the number of 
nodes that are sent to nodes of the same type. 

Proof. Let M : Vi — > V2 be any matching. Then 

w{M) = J2 w{v,M{v)) = H{v,M{v)) + ^Yl x{v,M{v)). 
veVi veVi veVi 

The first addend is a positive integer, while the second addend is strictly smaller 
than 1, because by Proposition l.(a) both A'^i and N2 have at most n — 1 hybrid 
nodes, and therefore ^.ygy x(^7 -^^('^')) ^ 2(n — 1). Therefore, Y2v<=v H{v,M{y)) 
is the integer part of w{M). This implies that w{M) ^ w{M') if, and only if, 

Yh{v,M{v))^ Y,H{v,M'{v)) 
and, if the latter are equal, also 

veVi veVi 
from where the statement clearly follows. □ 

Remark 5. If we restrict this alignment method to phylogenetic trees, the weight 
of a pair of nodes (^^1,^2) is simply \Cl{vi) A Cl{v2)\. This can be seen as an 
unnormalized version of the score used in TreeJuxtaposer [25]. 

Remark 6. Let N = {y,E) and N' = {V',E') be two tree-child phylogenetic net- 
works without out-degree 1 tree nodes. If they are isomorphic, the isomorphism 
between them is an optimal alignment of total weight 0. The converse implication is 
clearly false in general: a matching of total weight need not be an isomorphism. 
Consider for instance the optimal alignment between the phylogenetic trees de- 
scribed by the Newick strings (1, 2, 3) and (1, (2, 3)). 

But if there exists an alignment M (obviously optimal) between N and N' of 
total weight and if \V\ = \V'\, then M is a bijection between V and V that 
preserves the /x-vectors, because 

Y IfJ'iv) - lJ'{M{v))\ ^ w{M) = ^ iJ,{v) = li{M{v)) for every veV, 
and therefore lJi{N) = n^N'), which implies, by Theorem 1, = N' . 
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Given two 5-DAGs Ni = (^1,^1) and N2 = {V2,E2) together with their 
^-representations ^(iVi) and ijl{N2), the simple Algorithm 4 computes the total 
weight w{M) of an optimal alignment M between A'^i and in + |V2|)^) 

time. 

Algorithm 4. Given the ^-representations /i(iVi) and n{N2) of two S-DAGs Ni = 
(yi,Ei) and N2 = (y2,E2), compute the total weight of an optimal alignment be- 
tween Ni and N2. 

begin 

let G = {{Vi \ {Vi)l) U V2 \ (^2)l) E) be a complete bipartite graph with l^il + IV2I — 2n vertices 
for each xi eVi\ {Vi)l and X2 G V2 \ (V2)l do 

set weightlxi, X2] = abs{fi{xi) — fj.{x2)) 

if xi and X2 are not both tree nodes or both hybrid then 
add l/2n to weight[xi^X2] 
let M ■.Vi \ {Vi)l — > V2 \ (V2)z, be a minimum- weight bipartite matching of G 
extend M by sending each v £ to the leaf of with the same label 

return M 

return J2^^lweight[vi, M{vi)] 
end 

Example 9. Consider the tree T and the galled tree depicted in Fig. 3. The total 
weight of the matching between T and A^ that sends the root of T to the root of 
N and the node n of T to the tree node a of A is 

w{r, r) + w{u, a) + w{l, 1) + w{2, 2) + w{5, 3) = 0-hl-h0-h0-F0 = l, 

and hence, since /i(T) ^ ^(A), it is an optimal alignment: see Fig. 11. 




Fig. 11. An optimal alignment between the tree and the galled tree in Fig. 3. 



32 



Example 10. Consider the tree-child phylogenetic networks N and N' labeled in 
{1,...,5} given in Fig. 12. The /i-vectors of their internal nodes are given in 
Table 4. Table 5 gives the values of y) for every internal node x N and 
every internal node y of A^'. From this table, the optimal alignment marked in red 
in the table and depicted in Fig. 13 (where, to simplify the picture, the arrows 
joining each leaf of N to the homonymous leaf in A^' are omitted) is deduced: its 
total weight is 8. It is the only optimal alignment between these networks. 




Fig. 12. The networks N (left) and iV' (right) used in Example 10. 
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Table 4. /i-vectors of the internal nodes of the networks depicted in Fig. 12. 



A web tool that computes an optimal alignment of two tree-child phylogenetic 
networks with the same leaves and without out-degree 1 tree nodes is available at 
the Supplementary Material web page. 
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Table 5. Weights of the pairs of internal nodes of the networks depicted in Fig. 12. 



7 Conclusion 

In this paper we have presented two methods for comparing pairs of tree-child 
phylogenetic networks: a metric and an ahgnment algorithm. While the former 
quantifies the similarity between two networks in a way that allows to soundly es- 
tablish whether a network is more similar to a second one than to a third one, the 
latter allows the interactive visualization of the differences between two networks. 
They are respectively the first true distance and the first alignment method de- 
fined on a meaningful class of phylogenetic networks strictly extending the class of 
phylogenetic trees. Since the class of tree-child phylogenetic networks includes the 
galled trees, this distance and this alignment can be used to compare the latter. 

Tree-child phylogenetic networks have been recently proposed by S. J. Wilson 
as the class of networks where to look for meaningful phylogenies, but for the 
moment no reconstruction algorithm for tree-child phylogenetic networks has been 
developed. So, it remains an interesting open question to characterize the sets of 
sequences whose evolution can be explained by means of a tree-child network and 
to provide an algorithm to reconstruct this network, as well as to characterize the 
computational complexity of these problems. 

On the other hand, several reconstruction methods for time-consistent tree- 
sibling phylogenetic networks have been proposed by Nakhleh and collaborators. 
Since no true distance for these networks is known so far, it is an interesting open 
question whether our distance and alignment method can be extended to these 
networks or not. 
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Fig. 13. An optimal alignment between the networks in Fig. 12. 
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